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The microscopic approach to the description of the phase behaviour and critical 
phenomena in binary fluid mixtures is proposed. It is based on the method of 
£C) ' collective variables with a reference system. The physical nature of the order 

parameter in a binary mixture is discussed. The basic density measure (Ginsburg- 
ON ■ Landau-Wilson Hamiltonian) is obtained in the collective variable phase space 

which contains the variable connected with the order parameter of the system. 
It is shown that the problem can be reduced to the 3D Ising model in an external 
field. 
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1. Introduction 



X 

The study of phase transitions and critical phenomena in multicomponent fluid 
systems is very interesting from the theoretical as well as practical point of view. 
Whereas in one-component fluid systems only gas- liquid equilibria exist, three differ- 
ent types of two-phase equilibria have to be considered in fluid mixtures: gas-liquid, 
liquid-liquid and gas-gas equilibria. Despite the numerous experimental results now 
available [ |IJ- [ ||, theoretical achievements in understanding a microscopic mech- 
anism of phase behaviour and nonuniversal critical properties of multicomponent 
fluids are limited. Most of the theoretical studies devoted to this problem may be 
divided into three main groups: phenomenological theories, mean field approaches 



and integral equation methods. The phenomenological approaches [ §]- [ I 1J give 
predictions about critical exponents and scaling functions but no quantitative esti- 
mates of nonuniversal critical amplitudes are possible within this framework. The 
problem of the phase diagram sensitivity to the microscopic model also remains 
unsolved. The mean field theories [ [12[ and integral equation methods [ |13[- repro- 



duce different phase diagram types by varying the microscopic parameters but give 
qualitative estimates. 
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Of special interest are Refs.[ p2| - [ \Hj devoted to the study of both the universal 
and nonuniversal properties. They are based on the previously proposed approach 
to the study of the gas-liquid critical point in a one- component fluid., namely, the 
hierarchical reference theory (HRT) [ |28fl. On the microscopic Hamiltonian grounds, 
the HRT develops the renormalization group structure near a critical point. Recently 
this approach was also used for studying the 3D Ising model [ 2Jj. 

In spite of the doubtless success of the HRT the problem remains of constructing 
a theory that allows within a unified approach a complete description of the phase 
behaviour of multicomponent continuous systems beginning with the Hamiltonian 
and ending with the thermodynamic functions in the neighbourhood of the phase 
transition point. This program has already been accomplished in both the 3D Ising 



model [ |32]1 and a simple fluid near the gas-liquid critical point [ |33)-[|35|. Within the 



framework of the 4 model this approach has permitted to obtain the non-classical 
critical exponents and analytical expressions for thermodynamic functions . 

This theory has its origin in the approach based on the functional representation 
of a partition function in the collective variables (CV) space [ |3TJ|. First the method 



of CV was proposed for a study of the charged particles systemsf |30|-[ pT| and then 



it was applied to the second order phase transition theory [ |3"2"|. The point is that the 
statistical description of the phase transition is to be performed in the appropriate 
phase space specific for a certain physical model. Among the independent variables 
of this space there must be those connected with order parameters. This phase space 
forms a set of CV . Each of them is a mode of density fluctuations corresponding to 
the specific feature of the model under consideration. In particular, for a magnetic 
system the CV are variables connected with spin density fluctuation modes, for a 
one-component fluid - with particle density fluctuation modes. What is the content 
of the CV for a multicomponent system? We will answer this question below. 

Experiments have shown that fluids and fluid mixtures near the ordinary critical 
points belong to the universality class of Ising-like systems [ [|. Thus, a study of 
critical properties of multicomponent systems requires, on the one hand, an extension 
of the method worked out for one-component fluids and, on the other hand, their 



further development. In [ [36], |37| we developed the CV method with a reference 
system (RS) for the case of a grand canonical ensemble for a multicomponent system. 
Within the framework of this approach the phase diagram of a symmetrical mixture 
was examined in detail [ 3S|-[ 51]. Our previous study has been mainly restricted to 
the Gaussian model. In [ |37| we obtained the explicit form of the Ginzburg-Landau- 
Wilson Hamiltonian (0 4 model) for a binary symmetrical mixture in the vicinity of 
the gas-gas demixing critical point. In this paper we generalize the approach in the 
case of a non-symmetrical binary fluid system (the system of different size particles 
interacting via different potentials). 

The layout of the paper is as follows. We give a functional representation of a 
grand partition function of a two-component system in Section 2. The physical na- 
ture of the order parameter in a binary mixture is discussed in Section 3. In Section 
4 we construct the basic density measure (Ginzburg-Landau- Wilson Hamiltonian) 
with respect to the CV which include a variable corresponding to the order param- 
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ctcr. 



2. Functional representation of the grand partition function of a 
binary mixture 

Let us consider a classical two-component system of interacting particles con- 
sisting of N a particles of species a and iVj, particles of species b. The system is in 
volume V at temperature T. 

Let us assume that an interaction in the system has a pairwise additive character. 
The interaction potential between a 7 particle at r*j and a 5 particle at fj may be 
expressed as a sum of two terms: 

Ujs(rij) = Vv( r u) + (fiysirij), 

where Vv( r ) * s a potential of a short-range repulsion that will be chosen as an 
interaction between two hard spheres <r 77 and ass- </v( r ) * s an attractive part of 
the potential which dominates at large distances. An arbitrary positive function 
belonging to the L 2 class can be chosen as the potential 7 a(r). 

Further consideration of the problem is done in the extended phase space: in 
the phase space of the Cartesian coordinates of the particles and in the CV phase 
space. An interaction connected with repulsion (potential Vv( r )) * s considered in 
the space of the Cartesian coordinates of the particles. We call this two-component 
hard spheres system a reference system (RS). The interaction connected with an 
attraction (potential (p-ys( r ) ) is considered in the CV space. The phase space overflow 
is cancelled by introduction of the transition Jacobian. The contribution of the 
short-range forces to the long-range interaction screening is ensured by averaging 
this Jacobian over the RS. 

Then a grand partition function in the CV representation with a RS can be 
written as (for details see Appendix A): 

E = HrjHi, 

where S is the grand partition function of the RS. The thermodynamic and struc- 
tural properties of the RS are assumed to be known. Although it is known that 
mixtures with only repulsive interactions might undergo a phase transition [ 35fl, we 



assume that in the region of temperatures, concentrations and densities we are inter- 
ested in, thermodynamic functions of the RS remain analytic. Hi has the following 
form: 



— 1 



f P 
J (dp) (dc) exp [PpfPo + PPi c ° ~ ^7 Y,l V ( k )PkP-k + 

k 



W(k)c^ k + U(k)p % c_i\ J(p, c). (2.1) 

Here the following notations are introduced: 

P£ and are the CV connected with total density fluctuation modes and relative 
density (or concentration) fluctuation modes in the binary system. 
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Functions pf and \L\ have the form: 

= ^r(A*i+A*i), ^r = ^(A*i-A*5) (2-2) 



(the expressions for /ij are given in Appendix A) and are determined from the 
equations 

91,121 = W, (2.3) 



dpt4 



(N a ) - (N b ). (2.4) 



V{k) = (4> aa (k) + 4> bb (k) + 24> ab (k))/2 

W(k) = (4> aa (k) + Mk)-20 ab (k))/2 

U(k) = (<f>aa(k) - (f> bb (k))/2, (2.5) 

J(p, c) = J(duj) (c?7)exp i27T^2(u; k p k + 7 fc c fc ) J(u,i), (2.6) 



J(u, )7 ) = exp[EE^r :: E MW(0,...,0) 

n>li„>0 ' k 1 ...k„ 



x 



7fe 1 ---7^fe in+1 ---^J- (2-7) 

Index z n is used to indicate the number of variables 7g in the cumulant expansion 
(2.6). Cumulants are expressed as linear combinations of the partial cumulants 



M 7l ... 7n (see ( |5.3| )) and are presented for 71, . . . ,7^ = a, b and n < 4 in [ [37| (see 
Appendix B in [ g7j). 

Formulas ( |2.1| )-( |2~7| ) are the initial working formulas in our study of phase tran- 
sitions in binary fluids. 



3. The order parameter in a binary mixture 

A choice of the order parameter in multicomponent fluid mixtures is a serious 
problem because the character of the phase transition can be changing continuously 
from the pure gas-liquid transition to the mixing-demixing one. The question of the 
physical nature of the order parameter in binary fluid mixtures has been considered 
until recently from the point of view of both the phenomenological theory [ [TTJ, 
|Tl"f and the microscopic approach [ pH[] , [ |23fl , [ [57], [ [EZ| . Nowdays the commonly 
accepted idea is that both the gas-liquid and mixing-demixing phase transitions 
are accompanied by total density fluctuations as well as by relative density (or 
concentration) fluctuations. This is the only symmetrical mixture which exibits a 
complete distinction between these two processes [ [5jJ. However, most likely such an 
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"ideal" system does not occur in reality. In real mixtures the contribution from each 
type of the fluctuation processes changes along the critical curve. The evaluation of 
such contributions at each critical curve point is essential to the definition of the 
order parameter and to the understanding of the phase transition character in the 
mixture. It seems to us that in our approach the question of the physical nature of 
the order parameter has a consistent and clear solution. Here we shall briefly focus 
on it. 

Let us consider a Gaussian approximation of functional integral (|2.1| )-( f2~7| ). This 
approximation, also known as the random-phase approximation, yields the correct 
qualitative picture of the phenomena under consideration. As the result of the inte- 
gration over variables 7^ and u^, Hi can be rewritten as 

1 \dp) (dc) exp \po(l3f4 + Mi/A) + 



"1 



i_ a 

c (/^r + N2/A) - (M^Ki + M 1 (1) K 2 



1 



YJypkP- 



X 



where 



Au(k) + c^c_ % A 22 {k) + 2pzc_ n A 12 (k)] 



Ni = M [ 2 2) M^ - M£>M{ 



(1) lAV 



M 2 = M^M^ - M£>M{ 



(!) aAO) 



(3.1) 



An (A;) 
A 22 (k) 
A 12 (k) 



M. 



(2) 



A 



M. 



(o) 



M. 



(i) 



A 



A = M 2 w M 2 l 



r(°) n/f ( 2 ) 



(M«) 2 . 



(3.2) 



In order to determine the phase space of the CV connected with the order parameters 
we distinguish independent collective excitations by diagonalizing the square form 
in (|3.1| ) by means of the orthogonal transformation: 

A{k) m + B{k)^ k 



P% 



C(k)^ + D(k)^ k . 



(3.3) 



The explicit forms for coefficients A(k), B(k), C(k) and D(k) are given in Ap- 
pendix B. 



As a result, we have 



-n- 

k 



1 



A(*) 
DM 2 



(dri) {dQ exp [r ]0 {AM 1 + CM 2 ) + 
-(M 1 (0) K 1 + M 1 (1) K 2 )/(A(0))- 



2 T,( £ i^ k )VkV-k + £22(k)^_k) 
fc 



(3.4) 
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where 

e u (k) = -(A u (k) + A 22 (k) T yJ(A n (k) - A 22 (k)f + AA\ 2 {k)). (3.5) 



One of the quantities ( |3.5|) (or both) tends to zero as the critical temperature 
is approached. On the other hand, we have to find such a CV rj^ (or £r„) which 
is connected with the order parameter. Index k* must correspond to the point of 
minimum of one of the functions £\\{k) or e 22 (k) ( or both). These functions depend 
on temperature, attractive potentials 4>ys(k) and characteristics of the RS. The RS 
enters into fl3.5|) by cumulants M y s(k). M y s(k) can be expressed by the Fourier 
transforms of the direct correlation functions C^s{k) by means of the Ornstein- 
Zernike equations for a mixture. In [ [Q] the analytic solution of the Percus-Yevick 



equation for a hard sphere binary mixture was obtained. 




Figure 1. Coefficients £u{k) and e 22 (k) as functions of k for the NH3 — N 2 mixture. 
The solid curves represent the gas-gas demixing critical point (T = 413°K, rj = 
0.45, x = 0.5) and the dashed curves represent the gas-liquid critical point (T = 
373° K, rj = 0.12, x = 0.5) [0 

Coefficients £u(k) and e 22 (k) are studied both as wave vector functions at differ- 
ent values of temperature T, density rj and concentration x including the gas-liquid 
and mixing-demixing critical points (see Fig. 1) [ [RJ and as temperature functions 
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at k = (see Fig. 2) [ f£l|. The results show that branch Su{k) becomes a critical 
one no matter whether the system approaches the gas-liquid or gas-gas demixing 
critical point. Moreover, en(k) and £22 (k) have the minima at k = [40]. Hence we 
can draw the following conclusions: 

1. Branch £u{k) is always critical. 

2. Because Su{k) has the minimum at k — 0, the CV connected with the order 
parameter is 770 in the case of the gas-liquid critical point as well as in the case 
of the mixing-demixing phase transition. The particular form of rjo for each 
of these phenomena can be determined by means of the relations between the 
microscopic parameters, temperature, density and concentration of the system, 
e.g. by means of coefficients A, B, C and D. 

3. In the plane (p , c ) we have distinguished two directions: the direction of 
strong fluctuations 1] and the direction of weak fluctuations As a result, 
we can write the conditions for the binary mixture critical point in the form: 



d 2 n 



drj 2 

d 2 n 1 



0. 



drjod^o 
d 3 Q 



0. 



(3.6) 
(3.7) 
(3.8) 



where Q = — fcTlnS is a grand canonical potential. 
Now let us consider equations (|3.3|) at k = 0. From (|3.3|) it follows that 

rjo = ±D(0) Po tB(0)co 
£0 = tC(0) Po ± A(0)c o , 

where the upper sign corresponds to the case 

\A 12 \ = -A 12 (AD-BC=1) 

and the lower sign corresponds to 



\A 12 \ = A 12 (AD — BC — — 1). 



(3.9) 



(3.10) 



(3.11) 



On the other hand, ( |3.9| ) can be rewritten as 



r/ = ±p cos 9 ± c sin 9 
£0 = — Po sin 9 + Co cos 9. 



(3.12) 



Comparing ( |3.9j ) and ( |3.12| ) we can determine rotation angle 9 of axes 770 and £0 in 
the (po, Co) plane from the equation 



tan9 = — 
A 



(3.13) 
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r-0.4 




0.6 b'.'i i!6 '/'.'£ I'd ts 0.6 b'.'i 'i!b T.t 2.'6 '£s 



Figure 2. Coefficients en (A; = 0) and £22 (k = 0) as functions of the dimensionless 
temperature T at q = 1.0, a = 1.0, rj = 0.26 and x = 0.7 for different values of r 
(sa(k = 0) are obtained for a binary hard core Yukawa mixture [|4( 
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(in both cases ( |3.10|) and ( 3.11 )). In the case ( |3.11|) transformation (|3.12|) corresponds 
to both the mirror image with respect to the en axis and the rotation in the (po, en) 
plane. Thus, the proposed approach allows us, on microscopic grounds, to define the 
order parameter at each point along a critical curve and so to understand the phase 
transition character in the binary mixture. 




Figure 3. Density-concentration projection of the critical line of the model binary 
mixture at a = 1.0, q = 0.9 and r = 0.6(mean field approximation) 

Figures 3-5 show the (rj, x) projections of the (T, 77, x) critical surfaces of the 
model binary mixture for various values of microscopic parameters a, q and r. The 
arrows show the direction of the strong fluctuations (order parameter) along the 
critical curve in accordance with formula ( |3.13|) . Here the following notations are 



introduced: 77 is the packing density (77 = r\ a + r] b , rji = ■np i af i /Q 1 pi = (N)/V is the 
number density of species i), x is the concentration (x = (N b )/{N)), a = cr aa /a bb 
is the hard sphere ratio, an is the hard sphere diameter, q = —<f) bb (0)/ \ aa (O) | 
is the dimensionless "like" interaction strength and r = — 4> ab (0)/ \ aa (O) | is the 
"unlike" interaction strength (the form of <f>ij(r) is not specified). The critical surface 
is derived from the (f m f, V, x) surface by the equations for a binary mixture critical 
point (in terms of the Helmholtz free energy [ |I[). f m f is the Helmholtz free energy 
of a binary mixture in the mean field approximation (see Appendix C). 

The purpose of our further study is the calculation of the binary mixture be- 
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V 



Figure 4. Same as figure 3 at a = 0.9, q = 0.9 and r = 0.6 
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haviour in the vicinity of its critical points. Based on the Gaussian distribution 
(]3.1|) - (|3.2D we have determined the critical branch and, correspondingly, CV rjo con- 
nected with the order parameter. Now we shall construct the basic density measure 
with respect to CV r)z. As it is shown in [ |32|, in the vicinity of the critical point the 
basic density measure exists which includes higher powers of CV than the second 
power. We shall follow the program: (1) having passed from CV p^ and to CV rj^ 
and £g in ( |2.1| ), we shall integrate over variables £g with the Gaussian density mea- 
sure; (2) then we shall construct the basic density measure with respect to variables 
?7g (Ginzburg-Landau- Wilson Hamiltonian) . We shall restrict our consideration to 
the rf model. 




Figure 5. Same as figure 3 at a = 0.9, q = 0.9 and r = 0.8 



4. The microscopic Ginzburg-Landau-Wilson Hamiltonian in the 
vicinity of the binary mixture critical point 

We pass in ( |2.6| ) from CV p^ and to CV and £g: 



— i 



(drj) (df) exp [ry /ii + £0^2 - - Y.[ 7 lk r l-k P ( k ) + 

k 
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where 



/i 1 = /3(4u+ + CK), fa = P(Bt4 + Dfi) (4.2) 



P(k) = ^(A 2 V(k) + C 2 W(k) + 2ACU(k)) (4.3) 

R(k) = ^(B 2 V(k) + D 2 W(k) + 2BDU(k)) (4.4) 

Q(k) = ^(ABV(k) + CDW(k) + (AD + BC)U(k)) (4.5) 

J(V, = / {dx) m exp [z2tt E( W + £A) + E A»(x, 0) (4.6) 

k n > 1 



= { -^- E ^ n) (o.-)x 

n! 

■ ■ ■ #k n - in Xk n _ (in+l) ■ ■ ■ Xk n (4.7) 

Xk = M + C ^ H = B "% + D lk 

and are linear combinations of cumulants M^ n >. 

The square form in ( |4. 1|) is diagonal if 

Q = 0. (4.8) 



Taking into account formulas from Appendix B, equation ( |4.8j ) can be rewritten in 
the form: 

(1 - x)S aa - xS bb q - 1 



yjx(l ~ x)S ab 



0, (4.9) 



where Sij(k) is a two-particle partial structure factor of the RS. 
On the other hand, the square form in (|4.7|) is diagonal if 

= (4.10) 
and the last equality holds if ( |4.8|) holds. Equation ( |4.8| ) holds in the following cases: 



a symmetrical mixture: S aa = S bb , x — 1/2, q = 1 

a non-symmetrical mixture: (1 — x)S aa — xS bb = 0, q = 1 

a non-symmetrical mixture: q ^ 1. 
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Let us eliminate the linear term in ( |4.7| ) by the shift 

& = ££ + Mi (0) % 
and present D n (x, i?) as a sum of two terms: 

where 

D' n = D n {ti)+D n { X ,&), K = D n ( X ). 

Here D n {d) includes only the products of variables i}j:,D n (x, &) includes the mixed 
products of both variables and xti D n (x) includes only the products of variables 
X%- Let us consider the integral 



1 



J (dQ exp 
(-z2tt) 2 



k k 
1 



k 

where the following notations are introduced: 

M + = fa-R(P)M$ > \ A=J2K, 

n>3 



in ( 4.11 ) the prime on £j is omitted for clarity). 



If operator ^ is substituted for z27T$£, ( |4.12|) can be rewritten as 



(4.11) 



(4.12) 



T =Ef(l + {A)+i-(A 2 ) + ...), 



where 



(...) 



1 1 



E<; 2 11(^2 (*))"* exp(^) J (dfl exp [.M + £ 



. . . exp 



2 J ^ 



(0) 



Sf = exp(^) J] V^WMfW + lexp f - ( M+ ) 2M 2 0) 



V 2(i?(0)M 2 (0) + 1) 



Finally, after the integration in ([4.1|) over variables ^ we obtain 

S = ~ ~fAt J (dr])exp [fi^o - ^Y,VkV-k P ( k )] J (v), 

k 



(0)^ 



(4.13) 
(4.14) 

(4.15) 
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where A^ is the result of integrating over ^ which does not include x%- nas the 
form: 



r r (-i2n) 2 

J(V) = y(rfx)exp[z27r^ X ^-^7rM 1 (0)^xA- + ^^ M 2(0)X:XfeX-fc + 

k k k 

H27r) 3 „ . _ . (-i27r) 4 tr , n . 

'—^M 3 (0, . . .) XfeXfeXfe^+fe+fe + K —^M,(0, . . .) x 

E ^ ^XfeXfeXE,^^^! (4-16) 

kik2k$k4 

Here 

M n (0,...) = M(")(0,...) + AM n (0,...). 
Setting .M + = we get the following expressions for AM n (0, . . .): 

A? (1) 

AM i = 2(%E »(*') + ■■■ ( 4 - 17 ) 



M {2) (M {1 h 2 - 

AM 2 = ^p(k') + ^p(k'm\k-k'\) + 

^0 g sK*')S(*")sK|£ - k' - k"\) + . . . (4.18) 



AM 3 = ^MfM«^(A; 2 )^E^H^Mf x 

' E»( fc ')^(fe + fc3-fc'|) + ..- (4-19) 



where 



AM 4 = 3(M 3 (2) ) 2 ^(|4 + 4|) + 2Mi 1) Mi 3) ^x 

+ *3 + M) + ^(^i 2) ) 2 ^ X 

E^')ff(fe + fc4-fc'|) + ..., (4.20) 
k' 



l + R(k)Si 0) 
R(k)=R(k)(N), si 0) = M^/(N). 
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371. In this case 



An estimation of corrections AM n was carried out for the symmetrical mixture in [ 
M 2 = ! + M 3 = 0, M 4 = -2 + 5 2 , 



where 5{ are small values. 

Let us consider formula ( |4.3| ) for P{k). Substituting into ( |4.3| ) coefficients A, B,C 
and D from Appendix B, we get 



P(k) 



^{\(4>aa(k) + 4> bb (k)) + 



1 



AA\ 2 + (A u - A 



x 



22) 



lab(k)(A u - A 22 ) + (4> aa (k) - 4> bb (k))A 12 ] 



(4.21] 



where A n , A 22 and A 12 are functions of temperature. Using condition (|4.8| ), P{k) 
can be represented as 



P{k)V 



JL 

2V 



ka(k) + Mk) + 2 



A u - A 22 



1 



An - A 22 1 VI + k 2 
(4>ab(k) + K(4> aa (k) - 4> bb (k)) , 



(4.22) 



where 



q-1 (1 - x)S aa - xS, 



bb 



2Jx{l - x)S, 



ab 



P(k) takes both negative ( at small \k\) and positive ( at large \k\) values . In the 
region \k\ > B (see Fig. 6) we can integrate in ( 4.15|) over Xk anc ^ Vk w ith the 
Gaussian measure density as the basic one. As a result, we get in the approximation 
of the yf model 



~ Zf A«A r ' / d(rj) NB d{ X ) NB exp 



k 



-i2n) 



2! 



-M a (o)5> jEX _ 5 + 



(-i2vr) ; 
3T~ 



fcifc2fc3 

kik 2 k 3 k4 



k 

-Z27T) 4 



M 4 (0, . . .) x 



^1+^2+^3+^4 



(4.23) 



where in place of Mj(0, . . .) stand renormalized coefficients Mj(0, . . .): 

Mi(0, . . .) = Mi(0, . . .) + AMi(0, . . .) 



A "= II \/P{k)M 2 + l 

\k\>B 
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We can consider a set of k vectors, \k\ < B, as corresponding to the sites of a 
reciprocal lattice conjugated to a certain block lattice 77 with Nb block sites in 
volume V: 

{Nb) = v(Bhf = W(B fltt )'(s + n »(l-E» 

6tt 2 t] 

One may consider quantity B as the size of the first Brillouin zone of this block 
lattice. 

Next, two shifts are carried out in order to eliminate the cubic term in (|4.23f) [ 



34 



where 

A = -i2nM 3 

and 

- M 2 M 3 1(M 3 ) 3 
M, .I/, i— - + v A > 



M 4 3(M 4 ) 2 ' 
Then, (|4.23|) has the form: 



SoSf A« A" exp(£) j d{ri) N ° d(x') NB exp [hr,' - \ £ P(fc) 



x 

it 

, , (-il-Kf ~ , . _ , , (— i27r) 4 

fc ' fc 

M 4 (0,...) V x? Xs Xt xW- - - ~1 - ( 4 - 24 ) 

where 

M 2 (0) = M 2 (0) - ^^j^, M 4 (0, . . .) = M 4 (0, . . .), 

M 3 / - M 2 M 3 1(M 3 ) 3 \ - / 1~\ M1M3 

£ = — — Mi — + -t— ^7 IMj ft - -Mi + - 3 , 

M 4 V 2M 4 8 (M 3 ) 2 / 1 V 2 V M 4 ' 

ft = /ii - P(0)Mi + ^ 
After the integration over Xt m ( [4.24 ) we get 

S = C| exp^r?)]^)^, (4.25) 

where 

C = S Sf A«A"exp(^ + a N B )V2 NB '\ (4.26) 



a 4 



4! W 3 „ ,., 



X! ^ki^k 2 ^kz^ki^ki+k2+kz+ki, \ki\<B, (4.27) 
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do 



111 



1 /AW/4, 3 xi/4 /z 2 \ . 



N 
d 2 (k) 



LS41 
a 2 + P(k) 



d-2 



12 (N B ) 



\ \Sa 



Hz) 



a 4 = 36 



(N B ) 1 



(AT) 
/C(z) 



On 



15 



+ \m 



(N) 
K?{z 
--U(l,z)/U{0,z 

S 2 



3 (N) 



v / 2 n Mn/(A r }- 



(4.28) 
(4.29) 

(4.30) 

(4.31) 
(4.32) 

(4.33) 



Here U (a, z) is the parabolic cylinder function. 

As it is seen from ( f4.27| ), E±(rj) has the form analogous to the basic density 
measure of the 3D Ising model in an external field [ |32[] . But the main difference is 
the dependence of coefficients an, a 2 and 04 (see (|4.28|) - (|4.33|) ) on the microscopic 
parameters of the system. 



5. Conclusions 



We propose the microscopic approach to the study of phase transitions and 
critical phenomena in multicomponent mixtures. It is based on the CV method with 
RS. This method allows us to take into account the short-range and long-range 
interactions simultaneously. 

We consider the task of the definition of the order parameter in a binary mixture 
and show that it has a consistent and clear solution within the framework of our 
approach. 

After integration over CV £g (which do not include the variable connected with 
the order parameter) the basic density measure with respect to CV r]t (Ginsburg- 
Landau-Wilson Hamiltonian) is constructed. It is shown that the task can be re- 
duced to the 3D Ising model. 



Appendix A 

A grand partition function of a two-component fluid system in the CV represen- 
tation with the RS can be written as [ 3£| : 

E = E Ei, 



where 



00 00 b 

-0= e e n ex p 

Na=0N b =0 7=a 



AM 



(dr) exp 



7<5 ij 
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is the grand partition function of the RS; (5 — tj^ is the reciprocal temperature; 
(dT) = Ila,b^rAr 7 , <ir N-y = drldrl ■ ■ ■ df% is an element of the configurational space 
of the 7-th species; /i 7 is the chemical potential of the 7-th species in the RS. 

The part of the grand partition function which is defined in the CV phase space 
has the form of the functional integral: 

f P 
Si = j {dp)exp\PY,&PQn ~ ^77 E E Mk)p^p_z tS ]J(p a , t 5 ' 1 ) 

T 7-5 k 



Here, 



1) p[ is a part of the chemical potential of the 7-th species 

P 



Pi = Pj ~ Po + ^ E <f>ll( k ) 
k 



and is determined from the equation 

d In Ei 

dPul 



/i 7 is the full chemical potential of the 7-th species; 



(5.2) 



Pi 



■ipi is the collective variable of the 7-th species, the indices c and 



&i7 r k,y " r fc,7 

s denote the real part and the coefficient at the imaginary part of p^ 7 ; /o| and 

describe the value of /c-th fluctuation mode of the number of 7-th species particles. 
Each of p% and p% takes all the real values from —00 to +00. (dp) is a volume 
element of the CV phase space: 

(dp) = Yl d Po,-yU.' d P C k n d P%Y 

T MO 

The prime means that the product over fc is performed in the upper semispace; 
3) J(p a ,pb) is the transition Jacobian to the CV averaged on the RS: 



r b 

J(pa,Pb) = {du)Y[exp 



exp 



E 

n>l 



-i2?r) r 



x 



E E M 7l..-7n(4--- A>fc 1/n 

^-T"fcl...fc n 



&n -Tin 



(5.3) 



where variable z/^^ is conjugated to CV p^ 7 - M 7l „. 7n (fci, . . . , is the ra-th cumulant 
connected with S m ... yn (ki, . . . , fc n ), the n-particle partial structure factor of the RS, 
by means of the relation 



M. 



71-7) 



(/Ci, . . . , fc n ) 



N 7l . . . N ln S 7l ,„ Jn (ki, . . . , A^)^ +...+£„, 



where <% 1+ ... + £ n is the Kroneker symbol. 
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— * — * — * — * 

In general, the dependence of M 71 ... 7n (/ci, . . . , k n ) on wave vectors k\, . . . , k n is 
complicated. Hereafter we shall replace M 71 >7n (fci, • • • , fc n ) by their values in long- 
wave length limit M 71 ... 7n (0, . . . , 0); 

4) (f) 7 s(k) is the Fourier transform of attractive potential 7< $(r). Function (fi 7 s{k) 
satisfies the following requirements: <j)^s{k) is negative for small values of k and 
limu_> 00 <j>ts{k) = 0. The behaviour of (j) y s(r) in the region of the core r < a y s must 
be determined from the conditions of optimal separation of the interaction. For a 
very broad class of potentials the general form of (p-ys(k) is presented in figure 6. 




Figure 6. The Fourier transform of the attractive potential <f>^s{r) 



We pass in (|5.1]) to CV and eg (according to uj^ and 7g) by means of the 



orthogonal linear transformation: 

Pk = ~Y(Pk,a + Pk,b) ( 5 - 4 ) 

c k = -y(Pk, a - Pk,b)- ( 5 - 5 ) 

As a result, for Si we obtain formulas fl2.1p-(p77|). 
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Appendix B 

The coefficients A(k), B(k), C(k) and D(k) have the forms: 
A = v / 2|A 12 |[4A 2 2 + (A 11 -A 22 ) 2 -(A 11 -A 22 ) x 



\A xl -A 22 y + AA\ 2 ]-\ 
B = y/2\A 12 \[4A 2 12 + {A 11 -A 22 ) 2 + (A 11 -A 22 ) x 
^A u -A 22 y + 4A 2 2 ]-\ 

C = -^144 [^n - ^22 - J (A u - A 22 f + 4A 2 2 ] x 



2 (A 



12 



[4A 2 2 + (A n - A 22 ) 2 - (A n - A 22 )^(A U - A 22 f + 4A?,]- 1 , 
y/2\Ai 



D = -^^\[A ll -A 22 + ^(A 11 -A 22 )2 + 4A 2 2 \ x 

^ 1^12 J 



[4A 2 2 + (An - A 22 ) 2 + (An - A 22 ) J(A n - A 22 )2 + 4A 2 2 ]~ 1 . 



Appendix C 

The Helmholtz free energy of a binary mixture in the mean field approximation 
can be written as 

fmf fid fref fattri 

where fid is the free energy of a binary mixture of ideal gases, f re f is the free energy 
of a binary mixture of hard spheres [ [E| : 

f ref = F ref /(N)k B T = -1.5(1 - yi + y 2 + y 3 ) + (3y 2 + 2y 3 ){l - r])' 1 
+1.5(1 -y x -y 2 - y 3 /3)(l - vY 2 + (Vs - 1) m(l - V ), 

. 1 + a r] a a + r] b 
y 1 = A 12 ^^, y 2 = A12 = — , 

77 a v 



(1 — x)a 3 r] XT] 



Va = — 73 v— ^ 



x + (1 — x)a 3 ' x + (1 — x)a 3 

/attr = F attr / (N)kBT is the contribution due to attraction between the particles: 

1 n 

fattr = ~o 7fw —Tj T^T^ 1 ~ + 2X ( X ~~ + 

2 i *{x + (1 — xjcr) 

where T* = ksTa 3 | Q a(O) | _1 7r/6 is the dimensionless temperature. 
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